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ABSTRACT 


A  pertubation  expansion  is  formulated  for  the  three- 
dimensional,  nonlinear,  acoustic-wave  equation  with  dissipa¬ 
tive  term  describing  the  viscous  and  thermal  energy  losses 
encountered  in  a  cylindrical  cavity.  The  theoretical  results 
show  that  nonlinear  effects  in  sloshing  modes  are  strongly 
suppressed. 
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LIST  OF  SYMBOLS 


Q  phase  velocity  in  the  cavity 

dtf  =  cV  - 
Co  c^P A.?f  «t  P'P» 

^  effective  phase  speed  associated  with  a  standing  wave 

at  resonance 

Qj  specific  heat  at  constant  pressure 

-k  w/c’  propagation  constant  associated  with  a  standing 

wave 

peak  Mach  number  of  the  driven  standing  wave 
p  =  p  -  ~  acoustic  pressure 

p  P0  instantaneous  and  equilibrium  total  pressure  in  the 
field 

T<  absolute  temperature  in  kelvin 

particle  velocity 

particle  speed  of  the  T^h-order  perturbation  solution 
^  infinitesimal-amplitude  attenuation  constant 

oCnQo  temporal  decay  constant  of  a  resonance 

^  -('If  +  1)  /2  for  a  gas 

"Tf  =^P/Cvj  ratio  of  specific  heats  for  a  gas 

^ ^  instantaneous  and  equilibrium  densities 

velocity  potential 

yj  (angular)  grequency  at  which  the  cavity  is  driven 

]Mrl  (angular)  frequency  of  a  resonance 


W7/.V.. 
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I.  INTRODUCTION 

The  topic  of  finite  amplitude  acoustic  standing  waves  in 
sloshing  modes  of  a  cylindrical  cavity  is  interesting 
theoretically,  but  development  of  the  subject  has  not  been 
extensive.  The  purpose  of  this  research  is  to  present  the 
results  of  a  power  series  perturbation  approach  to  the 
problem. 


$ 

& 

& 

& 

& 


V  .• 
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II.  THE  NON-LINEAR  WAVE  EQUATION 

A.  GENERAL 

It  is  well  known  [l]  that  for  and  M«l,  where 
measures  the  fractional  loss  per  wavelength  and  M  is  the  peak 
Mach  number  of  the  source,  loss  terms  and  nonlinear  terms  in 
the  constitutive  equations  can  be  separately  approximated  with 
the  help  of  linear,  lossless  acoustic  relations.  The  force 
equation  appropriate  for  acoustical  processes  in  systems  for 
which  gravitational  effects  are  unimportant  is 


fjj-  +CU-7)U  +  |vp  =^U 


(2-1) 


Ky. 


where  (J  =  U  Li  is  the  particle  velocity,  ^is  the  instantaneous 
density  of  the  fluid,  P  is  the  instantaneous  total  pressure  in 
the  field,  and  the  operator  oC  symbolically  describes  those 
physical  processes  leading  to  absorption  and  dispersion.  We 
used  two  additional  equations.'  The  first  is  the  equation  of 
state  for  a  perfect  gas 


P  =  (>rT2 


(2-1-a) 


where  r  is  a  constant  whose  value  depends  on  the  particular 
gas  involved  and  Tk  is  the  absolute  temperature  in  Kelvin. 
The  second  is  the  continuity  equation 
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77  +  V-[C l+S)Uf  = 


=  o 


( 2-l-b) 


where  s  =  (  ^  -  ^  is  the  condensation  at  any  point  and  p 

is  the  equilibrium  density  of  the  fluid.  If  we  ignore 
rotational  effects,  then 


u  =Vl 


(2-2) 


where  1  is  the  velocity  potential.  Combination  of  Eqs.(2-1)- 
(2-2)  and  the  neglect  of  terms  of  orders  higher  then  M2 ,  MC°^/-fc)  , 
and  (  yields  a  quadratically  nonlinear  wave  equation, 


CC.‘[T+ -5: 


P  - 


u  * 


,  ,  Vf,  p  ,  u  ^ 

(cJj 


(2-3) 


where  Co2=(  )  (adiabatic),  p  =  acoustic  pressure  ,^A=(J  .  ^  , 

is  the  ratio  of  heat  capacities,  and 

cnx  =  cv*-4z 

The  left-hand  side  of  Eq.(2-3)  is  the  classical,  linear  wave 
equation  with  losses  pertinent  to  the  system  under  study. 
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The  right-hand  side  can  be  interpreted  as  a  forcing  function 
consisting  of  a  three-dimensional  spatial  distribution  of  phase- 
coherent  sources. 

In  a  second-order  perturbation  theory,  this  volume  forcing 
function  is  obtained  from  the  classical(first-order)  solution 
?1  of  the  acoustical  problem.  The  second-order  perturbation 
solution  P]_+P2  describes  the  non-linearities  ?2  resulting  from 
the  self  interaction  of  the  classical  solution  P-j_. 

Higher-order  perturbation  solutions  consider  the  interaction 
of  the  nonlinear  solution  with  itself,  and  the  forcing  function 
is  composed  of  products  of  both  classical  and  non-linearly 
generated  terms.  Thus,  if  a  system  is  driven  at  frequency  w, 
the  non-linear  term  in  Eq.(2-3)  will  force  the  existence  of 
all  integer  multiples  IflUf  of  the  driving  frequeney  and  the  full 
solution  must  contain  all  harmonics  of  the  input  frequency. 

In  a  closed  cavity,  each  of  those  nonlinearly  generated  waves 
whose  frequency  lies  near  the  resonance  frequency  of  a  standing 
wave  of  the  cavity  and  whose  assciated  spatial  function  matches 
that  of  the  standing  wave  can  be  strongly  excited  [2]  . 

As  far  as  the  author  has  been  able  to  determine,  there  has 
been  only  one  previous  study  of  this  system  published  in  the 
open  literature.  This  was  by  Maslen  and  Moore  in  1956  [3]  . 
Their  approach  resulted  in  a  series  expansion  which  did  not 
converge  if  the  relevant  nomal  mode  frequencies  were  integerally 
related.  Their  interpretation  of  the  quenching  of  the  nonlinear 
effect  was  based  on  the  "scattering  effect  of  the  wall". 
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We  feel  our  interpretation  based  upon  nonlinearally  generated 
volume  sources  stimulating  the  allowed  standing  waves  ot  the 
cavity  is  more  accessible  and  informative.  Further,  the 
mathematical  approach  developed  herein  appear  to  avoid  any 
difficulties  in  convergence.  Their  conclusion  that  high 
amplitude  monofrequency  transverse  oscillations  can  exist  is 
consistent  with  our  finding. 

B.  APPLICATION  TO  THE  CYLINDRICAL  CAVITY 

The  circular  cylinder  is  provided  with  a  "point"  source  of 
sound.  By  properly  positioning  the  cylinder  with  respect  to 
the  sound  source  it  is  possible  to  effectively  drive  the 
enclosed  air  into  various  modes  of  vibration.  The  rigid- 
walled  cylinder  forces  the  component  of  particle  velocity 
perpendicular  to  each  cavity  surface  to  vanish  at  the  surface. 
The  resulting  steady-state  solution  to  the  linear  wave  equation 
in  cylindrical  coordinates  is 

p  ^  f  q  ^  )  =  Xjfe  2  )  cc$  (7i<^ )  CCS  )  Jm C-^rmn 

'  '  7lin  (2-5) 


where  are  the  cylindrical  Bessel  functions  and  application 
of  the  boundary  condition  to  the  sides  yields 


where  a  is  the  radius  of  cylinder  and  are  the  arguments 

■tli 

of  the  exterema  of  the  "7!  Bessel  function. 

The  nomal  mode  frequency  is  dependent  on  and  ^vywn  , 


* 


fmr i 


+ 


-  G* 


(2-5-b) 


The  standing  wave  well  be  identified  by  the  ordered  integers 
(n,  m,  X  )  describing  its  spatial  dependence. 

1 .  Symmetric  Modes 

The  simplest  waves  in  cylindrical  coordinates  are  those 
that  depend  only  on  the  distance  r  from  Z-axis,  the  gradient 
takes  the  form 

where 

M  r 

0  (2-6-a: 


and  k=w/c  is  the  wave  number  or  propagation  constant. 
The  Laplacian  becomes 


v-7*  _  j )  * ,  a-  ,  i  JL  \ 

V  *  (  rf ‘  +  ^  yj ) 


( 2-6-b) 


The  suitable  trial  solution  appropriate  for  symmetric  modes  is 
the  linear  solutionCO,  m,  0)  of  frequency  w 


JL 

&c 


M  Jc  c-^)  tM* 


( 


The  velocity  potential  for  this  limit  can  be  defined  as 


(  % )  M  To  Cl)  cos  u>t 


.  u, , 

and  / r  is 

V*-'C 

-iL  =  ^  M  j/c^)a;s  wt=-jMJ,cc6wt  ( 

Because  we  assume  that  the  surfaces  of  the  cavity  are  rigid 
then  at  r=a  the  appropriate  boundary  condition  is  Jn'=0. 
Thus,  we  have  -&Q.  = j ^  . 

To  generate  the  second  order  solution,  we  first  note  that 

and 

MzJ*(l  +  cosz»t)  (: 


Next,  the  second  derivatives  of  Eqs.(2-1Q)  and  (2-11)  with 
respect  to  time  are 


jg  (■  %<j)  z  1  1  C{,S  Wt 


(2-12 ) 


3f  ^  &  ' 


=-c-?Vj>s^t 


(2-13) 


Substituting  Eqs.(2-2)  and  (2-13)  into  the  first  term  of  the 
right-hand  side  of  Eq.(2-3)  yields 


(2-14) 


and  the  second  term  of  right-hand  side  of  Eq.(2-3)  becomes 


c-  j:+ ■ +i*-  3i: )  ■ -( t-  JcV 

4-  ■* 


(2-15) 


since 


gVxJo'=  ( 3T*-J.*) 


(2-16) 


c.vjr-^ciT.'+ir-r) 


(2-17) 


With  the  use  of  Eqs.(2-14)  and  (2-15)  we  can  write  the  right- 
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hand  side  of  Eg. (2-3)  as 


s  { 


4 


Ml  h  imhi/j  +[-i 


■  COS  zurt  | 


(2-18) 


and  the  appropriate  inhomogeneous  wave  equation  for  the  second 
order  perturbation  solution  is 


2 .  Non -Symmetric  Modes 

The  non  planar  waves  in  cylindrical  coordinates  are 
those  that  depend  on  the  distance  r  from  Z  axis  and  the  angle 
from  X-axis.  The  gradient  takes  the  form 

i  2 


v  =  -£•  c  ) 


(2-20) 


where 


and  the  Laplacian  becomes 


v*  =  + 4^') 


(2-20-a) 


The  suitable  trial  solution  appropriate  to  the  forcing  function 
of  frequency  w  for  exciting  the  (n,  m,  0)  standing  wave  is 


where  Jn(^r)  is  the  Bessel  function  of  the  first  kind  and 
order  n.  The  velocity  potential  is  approximated  by 


(2-21) 


£=(-%-)  M  X  (•*’0  ccS  vut 


(2-22) 


and  U  */r  is 


iL  =  M(^  <*s  wt  -  §  X  ^ cc 6  )  (2-23) 


Thus , 


( )  =  M*X*  **  wt 


(2-24) 


(  JL  f= -i  mx  ■  -  -  wt)(x"+XlcfS 


( 2-25  ) 


and 


(-£)  =  £  rf(  I  XHt  ft  (  T,:+ J,H  )  -  JWX-iaS  Z^J 


(2-26) 


and  the  second  derivative  of  Eqs.(2-25)  and  (2-26)  with  respect 
to  time  is 


f  =  M  < i  ■ +**  - ^ ) coS 

di  n. 


(2-27) 


—  ( —  f 


_a£fM 

4  " 


r-Ki/x* 


+  [  X  -tf)  u 


vi _ i 


COS  Z'T) 


(2-28) 


Substituting  Eqs.(2-26)  and  (2-27)  into  the  first  term  of 


Eq.(2-3)  yields 


(*UT) 


(VJ^COS 


(2-29) 


(2-30) 


since 

J,'*  =  4  C  3»«  +  ) 

and 

^  Jo,*  =  X  (  3^,  ■+■  +  z  > 

a 

In  order  to  solve  the  second  term  of  Eq.(2-3),  we  can 

Co'V'j,2  =  ) 

G>x  cos  2,1*3  = 

Co  V  ^  =  (3thi  3m  J* 

“  Jtw-2  3i+*'3^ 


(2-31) 


(2-32  ) 


(2-33) 


(2-34) 


(2-35) 


(2-36) 


»  v*  •'*>  v-t  v  .•*  'sj  a  ■  t  *7*7n 

*.  v 

& 


m  r -  w ;  v?  w .  r?  -  .  r1.1  r.  ir,  jtj  y,  y.’  T,  IT.  -*.  y.  t*j  g\  ¥*.  r,  w;  y,  t^.  r.  J\r'^r'J 


With  Eqs . C2-30 ) ,  C2-38)  and  (2-39)  we  can  get  the  inhomogeneous 
wave  equation  for  the  second  order  perturbation  as 


|a 


f 


^  4  3t-a  )  "h  4  ( 


JL  -r2 - L  y  T_  )  CoS  aq1' 

^  J't  Z  J'TSc  / 


./  / 


f'f  COS  xwft  [-  7£  3*+*  +  x  vTrti 

4-  ^ 


g  Jinx  ^  3^r)^»”Xrw  X~  z 


X 


-  yiX]00^  J 


(2-40) 
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.  At-.V.-'. 


If  we  let  n=0  in  Eq(2-40),  it  reduces 
Substituting  n=l  into  Eq.(2-40)  yields 
forcing  term  resulting  from  the  (1,  m. 


to  Eq. ( 2-19 ) ,  as  it  must, 
the  equation  with 
0)  sloshing  mode. 


o)<»Z 


+ 


L^l  xilyt 

f  ’ 


-t- 


_L-r:l  /  JL _ L  \  ~yz  -j- — t" 

4.  Jo  v  71  i4  ^ If 


t 


<2  JsTi  zj  XX  C 


J- 

'  g 


(2-41) 


III.  METHOD  OF  SOLUTION 


Recall  that  the  equation  with  forcing  term  resulting  from 
theCl,  m,  0)  sloshing  mode  can  be  written  as  a  function  of 


frequency. 


I-TccS  XUft 


•  J?  +  X-X 


-  t-t-M 


Zj  16 


+{?XI-FJxJ<r(l:-T^i1C0S^ 


4*  d° 


(3-1) 


and  the  left-hand  side  of  Eq.(2-3)  is 


r  z3:/  ?  ,  J L  \  _  JL  -4-  -i-  f  -~ 

I  Co  C  >14*  M  )  Jjf1  '  PoCo 


11T  ^  ?t 


( 3-l-a) 


If  the  harmonics  of  the  frequency  at  which  the  cavity  is  driven 
are  not  close  to  any  resonant  frequencies,  we  can  ignore  the 
lossy  term.  So,  Eq. (3-l-a)  can  be  written  as 


r*»,Y— +  — —  )-  —  )-&_ 
Co*  3i') 


(3-2) 


a 


Let  us  assume  that  the  solution  of  Eq.(3-1)  can  be  expressed  as 


JL  =  ^ccs  «»t[u.c^)  +V»C1)  **5 


(3-3) 
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and  define 


nW  =  (-fer)C£?5  CCS*u/t  (3-4) 

0Ce* 

and 

,Bk-U-  =  M^LlriC'f^)  CXi^Zt/^t  (3-5) 

fyCc 

Combination  of  Eqs.C3-2)  and  (3-3)  yields 


A.  POWER  SERIES  METHOD 

In  order  to  solve  the  right-hand  side  of  Eq.(3-7),  we  can 
use  the  definition  [4J 


jvcf  = 


,  ))*ii  ^ 

(±1)  T- 

*  . 

(-[  y  p CV+yU+z&OirzrV _ 

f(V^)  -k! 


(3-8) 


where  P  C  ^  -r  I  )  =  J  ! 
1.  Coefficient 


a.  Substituting  V  -4L~  n+2  into  Eq.(3-8)  yields 


*  /  I  v^^V2- 

Tt ^ 

-fco 


4)/  -fe/ 


(  3-9) 


■  -£-=0 


2(1 >+&tx) 


( 3-10) 


where 


Bn*  (  0 


M  z-fef  Ar 


(  0! + •  z)  /  ]  (  Z'fH-feM)/  -fe/ 
<*  '  > 


(3-11) 
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b.  Substituting  V  =  n+2,,^  =  n-2  into  Eq.(3-8)  yields 


( &  vj/  x* ' 


(qf)4 -k-tz)  /  (m-tJd-Z')  !(x*‘k)l.-k! 


( 3-12  ) 


~Tr tz  ^  2Z  C*  (“*  ^  ) 


^+4) 


(3-13) 


where 


r  =  (  i  y*  c  ^  ) ! 

^  \  )  Crf+^-hZ,)  I  ( /rv+-£>  -z. 


C'l1+-fe4^)  1  (/TH-fenfc)  kt 


( 3-14) 


c.  Substituting  ^  =  n+l ,.((  =n-l  into  Eq.(3-8)  yields 


z_ 

(-i)^(^4^)  !  ^ 

(oH'-fei-i)/  Gu-HH )  /  (*nt£)? 


( 3-15) 


T™<i>T«hC1>  )** ( z  ^  I 


^m-k) 


(3-16) 


where 


D**  t 1 ) 


it 


(  2-nr)  -t-q-fe  ) 


("TH-fti-  I  )  (  )  •  C^#)!  "&• 


or 


where 


d.  Substituting  l)=44=n  into  Eq.(3-8)  yields 

-f  H  y-  (-oV-n-i-z-fe)! _ 


‘  *  ■"  (-hf  ^ 


1  <-1 )  =4_  UhS> 
-H 


_  /  .  v-fe _ (aon  )  / 


E*=(-') 


fj Z(Z7H^)!  ^! 


e.  Substituting  V =^( =n+l  into  Eq.(3-8)  yields 


X^K^l)  I 

-ft 


(xon~fe-t-*')  !  *fe! 


( 3-17) 


( 3-18) 


( 3-19) 


( 3-20) 


( 3-21) 
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where 


Substituting  ])=/(/,  =n-2  into  Eq.(3-8)  yields 


*feo 

(-if (;rm4-4)f  (i^f* 
[c 'jH-fe-*)/]  (jm*&-4)!  & 


( 3-27) 


(3-28) 


where 


0rw~(-i  ) 


4  (iD+ife-4)/ 


!  J  "(0-7H--4-4)  /  ■& ! 


(3-29) 


2.  Vn(y) 


With  the  help  of  £qs.(3-9)  through( 3-29 )  Eq.(3-7)  can 


be  written  as 


^  ^  _  (f)  j  + 1 J  V'ty 


=  jl.  CU-Ah  J(f  D' 

-|=o 


fiCntfc.) 


(3-30) 


-yov  •». 


where 


h  =  X  _  j_ 

n  z  $ 


C 3-30 -a) 


and 


a  = 


at 


( 3-30-b) 


Now,  let  us  assume  that 


0° 


1 


u/th4.) 


( 3-31) 


The  first  derivative  with  respect  to  y=-&r  can  be  expressed  as 


(3-32) 


<*1 


1 


and  the  second  derivative  with  respect  to  y-=^r  can  be  expressed 
as 


^  ^  '4=0 


( 3-33) 


Substituting  Eqs.(3-31)  throughC 3-33)  into  the  left-hand  side 


of  Eq.(3-30)  yields 
OO 


4=o 


X'Th-'k-l ) 
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Eq.(3-34)  must  be  equal  the  right-hand  side  of  Eq.(3-30).  Thus, 

( i)  1  X-fef 'A  +di  -  (gC# ' +  z  Ur'1'  A  E  THt) 

^  '  (3-35) 

where  k  =  0,  1,  2,  3  .  (3-36) 

If  Eq.(3-35)  is  substituted  into  Eq.(3-36),  then  Eq.(3-35)  can 
be  expressed  as 

(Xno  *  (rz&)  C»n-H  )  ^  1  )  (Xi\  lv  ^  Ej?o) 

a.-..  t(^)z(  ^+^)c  *■)  CU=-( tChi^D",  +  AE») 

dnx+(-^J  C^  +  5)(  ?  )  ft  La) 


-f  (zzj  (i-it+i+i  )  Cl*,-kn~  ~\i  -”fe+  x  Dmf*’  ^  Effe) 

(3-37) 


The  normal  component  of  particle  velocity  is  equal  to  zero  on 
the  boundaries.  From  the  application  of  the  boundary  condition 
at  the  rigid  boundary  at  r=a, 


dVr 

d  $ 


2*3  M'n'-Ne.Vi 


tr 


where 


Jt-H)  =  A  J/d.  (-&«■) 

Substituting  Eq.(3-36)  into  Eq.(3-38)  yields 


J-\'  \-2trl 


(3-38) 


+  . 


/aCnHh) 


=  O 


(  3-39) 


3. 


un(y> 

With  Eqs . ( 3-9 )  through ( 3-29) ,  Eq.(3-6)  can  be  expressed  as 


( 3-40) 


V\J 

where,  as  before,  A  =  ^/2-l/8  and  a  =  "~T~Z 


Now  let  us  assume  that 


=  k*  ( *1) 


2(th4-x) 


( 3-41) 


4 


■frc 


The  first  derivative  with  respect  to  y=-&T  can  be  expressed  as 

dll*  i  f  ,  „  I  ,'J (3-4  2) 

77  •=  x  /  b.«*  i 

d1 

and  the  second  derivative  with  respect  to  y=-j^rcan  be  expressed 

C<ntfe-z)(^t ir  ( i'jj 

~k=o 


as 


= 4  £  1 )  U  f 


( 3-43) 


Substituting  Eqs.(3-41)  throughC 3-4 3 )  into  the  left-hand  side  of 
Eq.(3-40)  yields 

Y_  ( -0 1  (')t+-fe)  -+ W-feH  )J  W  (^1) 

-feo 


-t-  L*.  ( 


zC^-z) 


( 3-44) 


Eq. ( 3-44 )  must  equal  the  right-hand  side  of  Eq.(3-40). 
Thus , 
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IvV 


V*  V. 


ft! 


I 


& 


(.^T ('T'+'t-^f  P'.-fe«  ■+■  -  (~ 7f  EE.u + f q--  vm 


^rHu.-fc-t  — )~A  E 


r.-ji- 


( 3-45) 


where  ~  1,  2,  3  . 

If  Eq . ( 3-45 )  is  substituted  into  Eq.(3-36),  then  Eq.(3-45)  can 
be  expressed  as 

ko+ck)1^-^'*'  =-7I& 


i  "no 


t6 

h  +(i)lc°i  -  if  k*  =~  i ©’w 

Lc  +  (£f(.'»T  U? 

:-±A  +^4  - 


TO 


I??*  4  (*t)  0+0  (6  St3^" f  Htu 


\n\  +  (JzJ  Q-n+zf  kt  _  it  &H+  f  E*i  AEir^*0 


UjL  +(^)X  OM  L^=-^G^+  +  H,>T^E^+K*fltk^. 


(3-46) 
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Application  of  the  boundary  condition  at  r=a  yields 


d  Ui 


d1tr 


=  4"  2_  (T;+i-^)  W  C  *  J  )  ~° 


where 


Substituting  Eq.(3-36)  into  Eq.(3-47)  yields 

+  (T!)U  (-Tt+oUj^j  ) 


.  0 


( 3-47) 


(3-48) 


IV.  CONCLUSIONS 

Recall  that  we  can  compute  CU  from  Eqs.(3-37)  through ( 3- 3 9 ) 
and  Eqs. (3-14),  (3-17),  (3-20).  Thus,  with  the  use  of  Eq.(3-31), 
we  can  get  V  .  From  Eqs.(3-46)  through( 3-48)  and  Eqs.(3-ll), 
(3-20),  (  3-2  3),  (3-26),  (  3-29  ),  we  can  compute  .  Thus,  with 
the  use  of  Eq.(3-14),  we  can  get  Un .  Therefore  we  can  compute 
Vn(y)  and  Un(y).  Let  us  calculate  the  finite  amplitude  effects 
resulting  from  the  nonlinear  distortion  of  a  forced  radial  mode. 
If  we  excite  a(0,  m,  0)  mode  and  obtain  the  pressure  at  the 
circumference,  then  -  jar*  •  Given  this  value  of  ,  use 

of  Eqs.(3-14),  (3-17),  (3-20)  give  the  quantities  C,  D,  E. 

From  Eqs.(3-35)  and  (3-36)  we  can  compute  (h'.fc  .  Thus,  we  can  get 
Vn(y)  from  Eq.(3-31).  Let  us  calculate  the  finite  amplitude 
effects  resulting  from  the  nonlinear  distortion  of  a  forced 
sloshing  mode.  If  we  excite  a(l,  m,  0)  mode  and  obtain  the  pre¬ 
ssure  at  the  circumference,  then  ^  =  .  Given  this  value 

of  ,  use  of  Eqs . ( 3-11) ,  (3-20),  (3-23),  (3-26),  (3-29)  give  the 
quantities  B.  E,  G,  H,  0.  From  Eqs. (3-45)  we  can  compute  . 
Thus,  we  can  get  Un(y)  from  Eq.(3-41).  Substituting  n=0  at 
radial  modes  into  Vn(y)  and  Un(y)  yields  Vg(y)=Ug(y).  Therefore, 
with  the  use  of  Eq.(3-3),  we  can  get  the  solution  of  Eq.(3-2). 

A.  RADIAL  MODES ( 0 ,  m,  0) 


J1  =  mx  (-&»*■) 
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=  zti  Voikos  xwc 


so  that 


=  C)m)  /  Jo  C  Jm,) 


B.  SLOSHING  MODES ( 1 ,  m,  0) 


=  1^|  31  (.'Mm  Ai* 


then 


/oC^+Vi^)^^]"5^ 


so  that 


liLl 

IM  ! 


)\m- 
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